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We present an ah initio study of the BFCO solid solution formed by multiferroics BiFeOs (BFO) 
and BiCoOs (BCO). We find that BFCO presents a strongly discontinuous morphotropic transition 
between BFO-like and BCO-like ferroelectric phases. Further, for all compositions such phases 
remain (meta)stable and retain well-differentiated properties. Our results thus suggest that an 
electric field can be used to switch between these structures, and show that such a switching involves 
large phase-change effects of various types, including piezoelectric, electric, and magnetoelectric ones. 
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Functional oxides attract attention because of their po- 
tential for designing materials tailored for specific appli- 
cations. A lot of work focuses on BiFeOa (BFO), one of 
the few compounds that is magnetoelectric (ME) multi- 
ferroic - i.e., displays coupled electric and magnetic or- 
ders - at room temperature [1 . Interest in BFO has been 
recently refueled by the discovery that an electric field S 
can be used to switch between two different ferroelectric 
(FE) phases of epitaxially-compressed films [2 . Such a S- 
switching has a number of functional effects associated to 
it, as the phases involved are markedly dissimilar in terms 
of cell shape (the switching thus implies a large piezo- 
electric effect) and magnetism (ME effect). Hence, BFO 
films offer the appealing possibility of obtaining phase- 
change functional responses of various kinds. Here we 
propose that some BFO-based solid solutions are ideally 
suited to this end, and present illustrative first-principles 
results for BiFei-ajCoa^Os. 

Materials- design aspects- The £^-switching in BFO 
films involves two phases [3j: one that is similar to the 
rhombohedral structure of bulk BFO and has a polar- 
ization P roughly along the [111] pseudo-cubic direction 
{R phase in the following); and a phase with a unit cell 
of very large aspect ratio {c/a ~ 1.25) and P roughly 
parallel to [001] {super-tetragonal or T phase). First- 
principles work has shown that these phases revert their 
relative stability as a function of epitaxial strain ^ |4j , 
and that there is a strain range in which both can ex- 
ist [4]. It has also been predicted that, even in absence of 
stabilizing fields, BFO presents many T phases that are 
local energy minima [5 . The theory is thus compatible 
with the observation that S fields can be used to switch 
between different FE phases of BFO. 

These results suggest that, to find materials in which 
the £^-switching is possible, one must look for compounds 
displaying a strongly discontinuous transition between 
two FE phases; further, the FE phases should be robustly 
stable and their polarizations point along markedly dif- 
ferent directions, so that it is easy to switch between 
them by applying properly oriented fields. An obvi- 
ous strategy is to look for chemical substitutions of 
BFO that may result in a morphotropic phase boundary 



(MPB) between the R phase of the pure compound and 
a second FE structure. Among many possibilities, the 
BiFei-a^Coa^Oa (BFCO) solid solution seems particularly 
promising. Note that bulk BiCoOa (BCO) is a ME mul- 
tiferroic that presents a super-tetragonal FE phase [6]; 
thus, BFCO is likely to display a R-T morphotropic tran- 
sition analogous to the one induced by epitaxial compres- 
sion in BFO films. 

The X-ray experiments of Azuma et al. [T confirmed 
that, as the Co content grows, BFCO moves from RtoT 
traversing a narrow region of presumably monoclinic (M) 
symmetry. (Similar results have been obtained for thin 
films [8 .) This is reminiscent of what occurs in prototype 
piezoelectric PbZri-a^Tia^Os (PZT), where the M phase 
found at the MPB [9 is characterized by its structural 
softness and large electromechanical responses [10 . How- 
ever, as far as we know, no enhanced response has been 
observed in BFCO, which questions the existence of a 
PZT-like M phase in this compound pT|. Note that this 
is encouraging in the present context: For £^-switching 
purposes, we would like BFCO's R and T phases to be 
relatively stiff (as opposed to soft) and stable. 

BFCO^s morphotropic transitions- We used the so- 
called "GGA+/7" approach to density functional theory 
(DFT) as implemented in the VASP package [12 , the 
calculation details being essentially identical as in pre- 
vious studies of similar materials [TJ. We worked with 
the 40-atom cell depicted in the insets of Fig. 1(b), which 
allows us to describe the R and T phases of interest [5] 
and vary the ratio of Co atoms x in steps of 1/8. Fig- 
ure 1(a) shows the formation energy of the phases in- 
vestigated as a function of composition; this is defined as 
= E — {l — x)E^YO —xEbco^ where E is the energy of 
a particular BFCO structure of composition x, and £^bfo 
and £^BCO are the ground-state energies of the pure com- 
pounds. Figures 1(b) and l(c-d) show, respectively, the 
results for the polarization (only for selected structures) 
and cell parameters. 

The six ponts at x = in Fig. 1(a) correspond to the 
stable BFO phases described in Ref. 5, whose proper- 
ties are summarized in Table [Tj The ground state is the 
well-known FE phase of the compound; we call it R-G, 
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TABLE I. Stable phases considered for BFO and ECO. We 
show the label for each phase (see text), space group (S. G.), 
polarization magnitude (/xC/cm^) and angle (relative to the 
perpendicular to the plane defined by the two shortest pseudo- 
cubic lattice vectors for the T phases, and to the pseudo-cubic 
[111] for the R phases), and energy (meV/f.u.) above that of 
the most stable phase. 



Material Phase 


S. G. 


Polarization 


BFO 


Ta-C 


Pc 


120 (23°) 


106 




Tc-C 


Cm 


150 (20°) 


103 




Tort-O 


Pna2i 


139 (0°) 


99 




Ta-C 


Cc 


145 (19°) 


96 




p-G 


Pnma 





27 




R-G 


R3c 


91 (0°) 





ECO 


Ra-G 


Pc 


82 (5°) 


52 




p-G 


Pnma 





51 




T-G 


PAmm 


167 (0°) 






noting its G-type anti-ferromagnetic (AFM) spin order 
(i.e., nearest-neighboring irons have anti-parallel spins). 
We also considered four FE T-phases with C-AFM or- 
der (i.e., parallel spins along the stretched lattice vector 
and anti-parallel in-plane). We label these phases ac- 
cording to their symmetry, which is not tetragonal for 
any of them: We have an orthorhombic phase (Tort), two 
monoclinic phases of type Ma (i.e., with P in the (110) 
plane [14]; we call them T4), and a monoclinic Mc phase 
(i.e., with P in the (100) plane; we call it Tc). Finally, 
we also considered the paraelectric phase p-G. 

To search for BCO phases, we substituted irons by 
cobalts in the six structures just described; then, for each 
of them we ran a short molecular dynamics (with ran- 
dom initial velocities to break all symmetries) followed 
by a full relaxation. We thus obtained the three solu- 
tions listed in Table I. All four T structures relaxed to 
the T-C phase known to be BCO's ground state [6 ; our 
computed structure (a = 6 = 3.70 A and c/a = 1.26) is in 
reasonable agreement with experiment {a = b = 3.729 A 
and c/a = 1.267 [6]) and previous DFT results [151 US]. 
The optimization starting from BFO's R-G phase led to 
a structure in which the R3c symmetry is slightly bro- 
ken to monoclinic Ma [faces of Oq octahedra lying on 
(111) planes form isosceles, instead of equilateral, tri- 
angles]; we call it Ra-G. Finally, BCO's p-G phase is 
analogous to BFO's. For intermediate compositions, we 
considered BFO's stable phases and studied all the in- 
equivalent Fe/Co arrangements at the perovskite B-sites. 
In all cases we ran a molecular dynamics followed by a 
full structural relaxation; the result of each such opti- 
mization renders a data point in Fig. [TJa). 

Many conclusions can be drawn from Fig. [l] Most 
importantly, we found that BFCO undergoes a strongly 
discontinuous transition between a R phase {c/a ~ 1) and 
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FIG. 1. (Color online.) (a) Formation energy vs composition 
for all studied phases. Labels as in Table [l] various instances 
of the same symbol at a given x correspond to different Fe/Co 
arrangements, (b) Polarization magnitude for the most stable 
structures (joined by thick line); a few others shown for com- 
parison; crosses correspond to rocksalt-ordered Ta and Ra 
structures. Insets: BFO (R-G) and BCO (T-C) structures, 
(c) Pseudo-cubic lattice constants (lower panel) and angles 
(upper panel) of the structures in panel (a); same-x config- 
urations given in order of increasing energy. Dotted vertical 
lines mark second-order R-R and T-T transitions (see text); 
solid line marks the first-order R-T transition. 



a T phase {c/a ^ 1.25) at x ^ 0.7. We also found that 
the R and T phases (as well as the p phase) are stable 
for all compositions [17 , a feature likely related to Bi's 
peculiar bonding properties [5] . Note also that the R and 
T phases retain their main features - i.e., values of po- 
larization and structural parameters - for all considered 
compositions and Fe/Co arrangements. Hence, accord- 
ing to our results, BFCO may allow for an £^-controlled 
switching between two distinct FE phases in a wide com- 
position range. Further, the robustly stable character 
observed for the R and T phases suggests there may be 
a relatively large experimental freedom to tune BFCO 
(e.g., by varying the composition or epitaxial conditions 
of a film) and optimize the switching. 

As regards the structure of the R phases. Figs, l(b-d) 
show no major changes occur when moving from BFO's 
R-G to BCO's Ra-G' Nevertheless, we identified a dis- 
tinct effect: For x > all the R structures present the dis- 
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tortion of the oxygen-octahedron faces described above 
for pure BCO. By inspecting the cases in which the Fe/Co 
arrangement is compatible with the rhombohedral 3-fold 
axis, we found this symmetry breaking renders a Ma 
structure, as in pure BCO. We thus label these phases as 
Ra-G in Fig. [l] where a dotted line at x ~ 0.07 marks a 
continuous transition between R-G and Ra-G. 

As for the T phases. Figs, l(c-d) show they evolve con- 
tinuously as the Co content increases, and become truly 
tetragonal (i.e., the pseudo-cubic lattice constants a and 
h become equal, and the angles turn 90°) for x ^ 0.8; this 
second-order transition to BCO's T-C phase is marked by 
a dotted line in Fig. 1. 

Our results thus confirm the R-T MPB reported by 
Azuma et al. [3; yet, the agreement between theory and 
experiment is far from perfect. At a quantitative level, 
our calculations place the MPB at x ~ 0.7, at variance 
with the value of x ~ 0.35 obtained by extrapolating the 
data of Ref. [7|to low temperatures. We think this devia- 
tion can be partly related to DFT's limitations to predict 
accurately the relative stablity of the R and T phases, as 
recently discussed for BFO [5]. The simulations, on the 
other hand, seem reliable when they predict that BFCO 
displays no PZT-like M phase acting as a bridge between 
the R and T structures at the MPB; rather, our results 
would be compatible with a R-T phase coexistence at the 
MPB region. Finally, more structural measurements are 
needed to confirm the monoclinic symmetry predicted for 
the R and T phases at intermediate compositions. 

BFCO^s phase-change properties - The R-T switching 
involves a large change in BFCO's unit cell, which results 
in a phase-change piezoelectric effect. It also involves a 
large change in polarization: We may switch between 
^ 53 /iC/cm^ {R phase) and ^ 165 /iC/cm^ (T 
phase), which may prove useful in the design of field- 
effect and other devices. Additionally, the R and T 
phases present very different dielectric and piezoelectric 
responses: For example, for x = 1/2 and a rocksalt Fe/Co 
order, we obtained an approximately diagonal static di- 
electric tensor for the Ra-G phase, with e^^^* ^46; in con- 
trast, the analogous Tc-C phase presents an anisotropic 
response with e^%** ^ 259, ejf^** ^ 122, and e^/^" ^ 16. (We 
observed similar trends for piezoelectricity. These results 
reflect a well-known fact for FE perovskite oxides: mak- 
ing P rotate is energetically less costly than changing its 
magnitude [10 .) Hence, the R-T switching also allows 
for a large dielectric and piezoelectric tunability ^18^ . 

The R-T switching also involves a change in the spin 
order, which moves from G-AFM to C-AFM. If any, the 
net magnetization of such AFM structures will be a small 
one arising from spin canting (see Ref. |4] for representa- 
tive results for BFO); thus, the associated phase-change 
effect will be tiny. On the other hand, work on BFO OH] 
shows that the magnitude of the exchange interactions 
varies considerably between the R and T phases. We 
found that such a differentiated behavior also occurs in 
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FIG. 2. (Color online.) (a) and (b): Energy differences be- 
tween spin orders for representative BFCO structures (see 
text); R phases in (a) and T phases in (b). Circles show 
DFT results; crosses correspond to model Hamiltonians in- 
cluding the J couplings defined in insets (see text), (c) Tem- 
perature dependence of AFM order parameter (G-type for 
R and C-type for T). (d) Temperature dependence of the 
total magnetic susceptibility (solid symbols), and of its com- 
ponents [parallel (||) and perpendicular (±); open symbols] 
with respect to the direction of the AFM order parameter. 



BFCO, resulting in markedly different Neel temperatures 
(Tn's) and response properties. 

To investigate BFCO's magnetic properties, we con- 
sidered the most stable R and T configurations at x = 0, 
1/8, 7/8, and 1, as well as the lowest-energy structures 
with a rocksalt Fe/Co arrangement at x = 1/2 [19]. We 
computed the energies of the following spin orders: fer- 
romagnetic (FM), C-AFM, G-AFM, and A-AFM (as ob- 
tained from G-AFM by making spins parallel in-plane); 
the results are shown in Figs. 2(a) and 2(b). Most no- 
tably, we found that the structure type determines the 
hierarchy of spin arrangements: Independently of the Co 
content, all the R phases have the G-AFM order as the 
most stable one, followed by C-AFM, A-AFM, and FM; 
in turn, C-AFM is the ground state of all the T phases, 
followed by G-AFM, A-AFM, and FM. From a model 
Hamiltonian perspective, this implies that the magnetic 
couplings will depend strongly on the atomic structure 
{R or T) and weakly on the chemical details (Fe/Co ra- 
tio and spatial arrangement). Indeed, we were able to 
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capture the essence of our DFT results in the two simple 
Heisenberg models, one for all the R phases and one for 
all the T phases, described in the following. 

We write the energy of a spin configuration as £^ = 
-\- 1/2 Jij Si • Sj^ where is a reference energy, 
Jij is the exchange coupling between Fe/Co atoms i and 
j, and Si is the three-dimensional spin of atom i. For 
simplicity, here we assume \Si\ = 1; the reported Jij 
values have been calculated accordingly. 

In the R phases, any Fe/Co atom has six Fe/Co first 
nearest neighbors (fnn's) that are roughly equivalent 
(they are equivalent if the symmetry is exactly rhom- 
bohedral). Thus, we can tentatively take Jij = Ji as our 
single fnn interaction, and also neglect further couplings. 
From our results for the pure compounds, we obtained 
Jf^ = 38.0 meV and Jf^ = 38.1 meV, and from the 
results at X = 1/2 we got Jf^ =38.1 meV. The minimal 
Hamiltonian thus defined predicts the energies shown as 
crosses in Fig. 2(a); it clearly captures the essential mag- 
netic interactions in BFCO's R phases. 

For the T phases, we take advantage of the (approx- 
imate) tetragonal symmetry and consider only two fnn 
couplings: in-plane Ji and out-of-plane Jy . Note that 
by choosing Ji > and Ji' < we would correctly repro- 
duce the C-AFM ground state, but wrongly predict the 
A-AFM order as the least favorable one. Hence, we need 
to introduce a new interaction J2 between second nearest 
neighbors (snn) [see inset in Fig. 2(b)]. With this model 
we obtained the results shown as crosses in Fig. 2(b); 
the computed parameters are (in meV): J^^ = 30.9, 
Jf^ = 35.0, Jf^ = 40.8, J^;^ = 2.3, JfP = 3.3, 
Jf,^ = 4.5, J^^ = 2.3, Ji^ = 2.0, and J^^ = 1.6. Note 
that all the computed Ji/ and J2 constants are positive, 
implying that all the out-of-plane interations are anti- 
ferromagnetic in the T phases of BFCO. Ultimately, the 
J2 couplings prevail and render the out-of-plane parallel 
spin alignment of the C-AFM ground state. 

We solved these Hamiltonians using standard Monte 
Carlo techniques. Figures 2(c) and 2(d) show represen- 
tative results obtained for x = 1/2 systems with Fe/Co 
atoms randomly distributed in a 10x10x10 periodically- 
repeated simulation box. We found that the R and T 
cases present markedly different Tn's of about 675 K and 
425 K, respectively. These results are compatible with 
the experimental values for BFO (643 K) [1 and BCO 
(470 K) 0, and refiect the weak out-of-plane interac- 
tions that result in a lower Tn in the T case. Figure 2(d) 
shows the magnetic susceptibility: At room temperature, 
the total susceptibility of the T phase almost doubles the 
result for the R phase; hence, we predict that the R-T 
switching will have a clear magnetic signature. 

We have thus shown that the BiFei-^jCoa^Os solid solu- 
tions are likely to present large electric-field-driven phase- 
change effects of various types (piezoelectric, electric, 
magnetoelectric) . We hope our predictions will attract 
interest towards such promising materials and effects. 
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